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Abstract—The main objective of this paper is to take a 
look at Electromagnetic Field Transformation from the 
Wave Equation point of view, which is slightly different 
from the traditional one based on the Transformation 
Optics (TO) approach, as we explain below. We begin 
with a review of Transformation Optics (TO), which 
provides a novel paradigm for designing a whole host of 
electromagnetic devices by using geometry 
transformation. The underlying concepts of TO as well as 
its theoretical foundations are well established, thanks to 
the pioneering works of a number of leading researchers, 
which date back to 2006. In this paper we examine the 
realizability issue of the materials, prescribed by the TO, 
from a practical point of view and identify some 
difficulties we may encounter in the process of fabricating 
electromagnetic devices such as cloaks and microwave 
antennas following the TO prescription. Next, we 
introduce a Wave-equation-based (WEB) approach to 
field transformation to examine whether we could 
mitigate the problems encountered with the TO that we 
mentioned above. Finally, we introduce the Generalized 
Scattering Matrix (GSM) concepts of Field 
Transformation, and we use it to modify the performance 
metrics of the devices in order to render them realizable 
by using available materials, without the need to use 
Metamaterials that are often lossy, narrowband, and 
polarization-dependent. 


Index Terms—Electromagnetic Field Transformation; 
Transformation Optics (TO); Wave Equation; Scattering 
Matrix Formalisms 


I. INTRODUCTION 


Perhaps no other topic has caught the recent 
attention of physicists and microwave engineers 
than the subject of Transformation Optics (TO) 
aka Transformation Electromagnetics (TEM), as 
the number of publications in the scientific 


literature [1-25] would indicate. This is because 
the TO provides an elegant theoretical approach 
for designing a variety of electromagnetic devices, 
e.g., Cloaks, flat lenses, reflectarrays, to name just 
a few, based on a fresh new way of looking at the 
problem that relates the electromagnetic fields in 
one coordinate system to another. For instance if 
we wish to design a flat lens following the TO 
algorithm, we can start with a plano-convex lens 
in free space and transform it into a rectangular 
shape with material parameters dictated by the TO, 
which maps the electromagnetic fields from one 
space to the other so that the two sets of fields are 
equal, and the performance of the flat lens exactly 
mimics that of the convex lens, provided we can 
realize the requisite material parameters of the flat 
lens, and of the medium surrounding the lens, in 
accordance with the dictates of the TO algorithm. 
Since the TO paradigm is well researched, and the 
algorithm itself is based on solid theoretical 
foundations of co-ordinate transformation of 
Maxwell’s equations [1, 2], it is legitimate to 
wonder if we really need another contribution on 
the subject, unless it adds something new and not 
just rehashes what is well known already. Our 
objective in this paper is indeed to introduce a 
different way of looking at the TO approach, 
focusing more on the practical aspect of realizing 
the electromagnetic devices whose geometries we 
choose, and whose performance we desire to 
emulate such that it corresponds to that of an 
existing device with a different geometrical shape. 
We have underscored the word “performance” 
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because, unlike TO, we may choose to work with 
certain selected “performance metrics” of the 
devices rather than attempting to render the 
Electric and Magnetic field distributions of the 
original and transformed devices to be strictly 
identical, as we do when we follow the TO 
algorithm. To follow our chosen path, we also 
deviate from the TO strategy, which works with 
Maxwell’s equations for both the Electric and 
Magnetic fields, and transforms these field 
equations from one space to another by using the 
rules of coordinate transformation to determine 
the material parameters of the geometrically- 
transformed device. Instead, we find that it is 
expedient for us to work with the “wave equation” 
satisfied by either the Electric or the Magnetic 
field, though we prefer the E-field wave equation 
because it enables us to work with the 
permittivity parameter € , and leave the 
permeability u unchanged. Working with e€ alone, 
if possible, would be significantly advantageous 
from a practical point of view, since realizing 
magnetic materials typically called for by the TO 
[3] has remained an elusive problem, even when 
the material is isotropic (anisotropic simply 
exacerbates the problem). Turning to 
Metamaterials (MTMs) does not really circumvent 
the realization problem, since the required MTMs 
can be lossy, narrowband, dispersive and highly 
polarization-dependent [43-48]. 


Following the discussion of the wave-equation- 
based (WEB) approaches, one of which is based 
on conformal mapping via the Schwarz- 
Christoffel transformation applied to the wave 
equation, we will turn to a slightly different 
formulation of the problem based on the 
Generalized Scattering Matrix (GSM) approach 
[26], and show how we can choose the 
performance metric of the electromagnetic device 
we wish to design to be slightly different from that 
of the TO, so that we can work with either 
available or with conveniently realizable materials, 
as opposed to MTMs typically called for by the 
TO. We will include a number of practical 
examples to illustrate some practical applications 
of the techniques we describe in this paper, 
reminding the reader that the practical aspects are 
indeed the main focus of our work. 


II. Review of TO and Introduction of 
Wave-equation-based (WEB) 
approach 


In this section we present the wave-equation- 
based (WEB) approach to field transformation, 
where our present goal—to be modified later in 
Sec.4 and explained why—is to preserve the E 
and H fields under coordinate transformation. 
Traditionally, TO accomplishes the field mapping 
by using the transformation of Maxwell’s 
equations from unprimed to prime coordinate 
systems as follows. 

VX E' = —jwu' H’ (la) 

VX H = +jwe'E' (1b) 

and shows that if we use the material parameters 
given by 
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then the E- and H-fields map one-to-one from one 
coordinate system to another, and we can achieve 
the same performance from the transformed 
device as we did form the original one, whose 
geometry we transformed from the original to a 
more desirable one. Experience has shown that the 
realization of the TO-dictated materials, examples 
of which are shown in Fig.l for the relatively 
simple geometry of a cylindrical cloak, is a 
daunting task, which is yet to be tackled 
satisfactorily even for such a simple geometry, let 
alone for a complex structure such as an airplane 
or a missile. 


Before we discuss the WEB approach, we turn to 
examine another fundamental limitation posed by 
the TO paradigm when we attempt to reduce the 
thickness of the cloak t = R, — R, to realistic 
values, e.g., a small fraction of the operating 
wavelength. We note that the material values 
(€, 4) do not change near the outer boundary of 
the cloak, located at p = R3, regardless of 
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Fig.1. (a) Schematics for a cloaked PEC cylinder with 
R» = 2R,; (b) material parameters for an all-angle, all- 
polarization cloak; (c) material parameters for a 
cloak—2D TM-polarization and normal incidence case. 


whether the thickness t of the cloak is large or 
small; and,for that matter, neither do the 
behaviors of (e,u) of the cloak in the 
neighborhood of p = R4 , where the above 
parameters approach extreme values. While we 
can attempt to partially mitigate this problem by 
imposing a chosen limit on the values of these 
parameters, we cannot change the fact that the 
relative (e(p), u(p)) must reduce from extreme 
values at the surface of the cylinder to moderate 
ones within a relatively short distance, which 
poses realizability problems in practice. Practical 
realizability of thin invisibility cloaks have not 
met with too much success in the past for this 
reason, and it is unlikely that the thickness issue 





Fig.2. Field distributions for a PEC cylinder: (a) a TO 
cloak with ideal material parameters; (b) a thin TO 
cloak with corresponding ideal material parameters; (c) 
a 3-layered TO cloak in which the medium parameters 


at each layer correspond to those of the ideal thin cloak 
in (b). 


will be resolved anytime soon if we continue to 
impose the above invisibility criterion on the TO- 
based cloak designs. Fig.2 shows the deterioration 
of the performance of the cloak (see Fig.2(c)) 
when we design a thin three-layer cloak by using 
discretized values of the continuous cloak 
parameters. It becomes evident, when we compare 
the performance of the fabricated thin cloak with 
the ideal one, that the wavefront of the field 
becomes considerably distorted when we use a 
thin cloak. We will address this problem later 
when we introduce a different performance metric 
for the cloak (or blanket) in Sec.3, and discuss a 
way we can use a thin coating to reduce the 
scattering form a target, albeit in the 
backscattering hemisphere. 


For the WEB approaches, which we present 
below, and which deal with coordinate 
transformation of the wave equation, we will 
initially limit ourselves to two-dimensional (2D) 
problems, deferring to examine the more complex 
3D case until later. We can gain much insight 
from the study of the 2D cases and draw some 
conclusions that are also valid and applicable to 
the three-dimensional case. We hasten to point out 
that similar two-dimensional cases have also been 
studied extensively by other workers, both 
theoretically and experimentally, and the 
cylindrical and carpet cloaks are but two examples 
of the 2D problem that have been researched 
earlier [3, 4]. 


Taking a cue from the literature we investigate 
the cases of conventional and carpet cloak designs, 
for an infinite cylinder and a triangular bump in an 
infinite PEC plane, respectively. We assume, also 
in common with the existing literature pertaining 


to these geometries, that the fields are independent 
of the longitudinal direction z, and that all the 
fields can be derived from the F,-component of 
the Electric field, which is the only non-zero 
component of the Electric field, from Maxwell’s 
equations 
V x E(x, y) = —jopH (x, y) 
(4a) 
Vx A(x, y) = +jweE (x,y) 

(4b) 
by setting E = ZE,. This case has been referred to 
in the physics literature [5] as the TE, case, 
although in the conventional waveguide 
terminology [27, 28] this case is termed TM, 
instead, and we will follow this latter terminology 
in this work. 


At this point we set u = My at the outset, and we 
will soon see how, in contrast to the TO paradigm, 
the WEB approach enables us to do this in a 
straightforward way. As we well know, finding 
materials with u values called for by the TO has 
been a severe roadblock and we wish to avoid this 
problem altogether, or to rely on MTMs for their 
realization. We wish to see what we can do 
instead with €-only materials for the cloaks and 
study how the WEB approach provides a very 
simple way to do this. Furthermore, we assume 
that € is isotropic (for 2D problems), and that it is 
invariant of the longitudinal coordinate z. If 
successful, this would obviously facilitate the 
realization of the cloak, since anisotropic 
materials with arbitrary properties, typically 
dictated by the TO, are very difficult to realize in 
practice. 

Next, we return the two Maxwell’s curl 
equations and eliminate H, by using 


1 
H(x,y) = in E(x,y) (5) 
0 
and substituting it into (4b), to get 
V x (V x E(x, y)) — w’ uoEE(x, y) = 0 (6) 


Using vector identities, we rewrite (6) as 
V(V : E(x, y)) — V?E(x, y) — w7 Up €E(x, y) = 0 


(7) 
From Maxwell's Equations, in a source-free 
region, we also have 


V: D(x, y) = V: [eE(x, y)] = Ve: E(x, y) + 
EV:E(x,y)=eEV:E@y)=0_ (8) 





24 











2b 2b 


Fig.3. The original (p,q,z) system (left) and the 
transformed (p’, p, Z) system (right) coordinate systems. 


which enables us to simplify (7) to: 
(V7 +k*)¥ =0 

(9) 
which is the Helmholtz equation (wave equation) 
satisfied by ¥ = E,. We wish to see how this 
equation transforms from one coordinate system 
to another, so that the large cylinder covered by an 
appropriately designed cloak behaves similar to 
the small one surrounded by free space instead. 


Thus we begin with the problem of deigning 
a cloak to reduce the scattering from a circular 
cylinder using the WEB approach by looking at a 
transformation of the two coordinate systems that 
relates the first system, namely the cylindrical 
region of 0 < p < b to the second, which is an 
annular region a < p’ < b, where p and p' are the 
radial coordinates of the two systems, respectively; 
ais the radius of the cylinder; and b is the outer 
radius of the cloak region, as shown in Fig.3 (also 
see [3], where R, and R, correspond to a and b, 
respectively). A linear transformation that relates 
the two coordinates is 


p'="“pta: 


b' =; z' =z (10) 
where @ and Z, respectively, are the angular and 
vertical coordinates in the original system (which 
is inherently free space) and @’ and z’ are the 
corresponding angular and vertical coordinates in 
the transformed system, respectively, which 


includes the cylinder and the cloak. 


As mentioned earlier, we wish to look at the 
problem from the perspective of wave equations, 


more specifically, the Helmholtz equation 
satisfied by the E-field: 

V7E,+k6E, = 0 (11) 

For the cylindrical system under 


consideration here, we rewrite (11) as 


1 3 dE,(p,o) , 1 0°EZ(p,o) 2 B 
pop 0p + p? agp? + kõEz(p, p) =0 
(12) 


or simply, 


10 0E,(Pp, p) 1 0°E,(9, p) 
pap’ ap p> dg? 
W*Up€oEz(p,¢) =9 (13) 
In the transformed coordinate p’, Q, Z system, 


the Helmholtz equation of the E-field can be 
derived to read as follows 


1 Cm joe) z(p')) 
(p' —a) are dp' 


1  07E,(p'.p) , ko? 
(p'-a)> 0? = 


b 











zE,(p',6) = 0 (14) 


and since we have chosen the permeability u 
to be that of free space, i.e., Ug (14) can be written 
in the form 


1 TE q) EN H) z(p P)) 
(p' —a) are dp' 





1 0°Ez(p',h) , @* Ho€o 
(p'—a)? Ap? Te o (es z E zp’, op) = T 0 (15) 





where the cylinder of radius a is mapped to a 
cylinder of limiting radius (tending to zero). In the 
general case, we denote c to be the radius of the 
small cylinder as shown in Fig.4, which yields 


p' =o; 


For this scenario the Helmholtz equation 
satisfied by the E-field can be rewritten as 


F z(p', ee 


) _ bra ac, 
P =P tb a 


= z =Z (16) 
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1 ð , 
(p -or (r - 
b-c 
1 0°E,(p’, 
—— A + aL Fe ($) = 0 (17) 
(p'—b5-5) (=) 
Next we shift the origin of the p’, ġ, z system 
by defining 


=o eee 
=p -b7 (18) 








Fig.4. A general form of the original (p, p, Zz) system 
(left) and the transformed (p’,@,z) system (right) 
coordinate systems. 


and rewrite (17) as 


a—c 
1 ð dE, (p, +b e) P 
pı 9/1 a ð pı 


1 0°E, (pa F bi, p) 


p1’ ap? 


p ai. (m +b 


=) 


In the p’ or p4 system, the medium 
parameter are UM, and €,€ 9, where €, is yet 
to be determined. Let 





=.) = 0(19) 








k2 
kf = kĝer = -a (20) 
(Fe) 
which yields 
S 
er = (Ba) 21) 
Then E, satisfies: 
a—c 
1 ð JE, (pi + bp —, P) 
pı 9/1 a Op, 
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5 1 0°E,(p, + br=) 
pr? ao: o 
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which is the form of the wave equation in the 
pı system. The choice of the material parameters, 
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(c) Relative permittivity 


Fig.5. Mesh distribution in: (a) original system in free 
space; (b) transformed system; (c) permittivity 
distribution in transformed system derived by using 
locally linear approximation of the transformation. 


given in (21), renders the form of the two wave 
equations (13) and (22) in these systems the same, 
and enables us to achieve our goal. 


For the limiting case c > 0, we need to set 


e= (>) Eo (23) 


which is the case of the ideal cloak with no 
scattering. For reduced scattering we choose a 
non-zero but small c. 


If we choose the outer radius b of the cloak to be 
2a, and if we let c = 0, then we see from (23) that 


the €, value in the p’, ġ, z system would be €, = 4. 


It is interesting to note, by referring Fig.1(c), that 
this is also the material parameter for €, derived 
by using the TO and presented in [5], where fg 1s 
seen to equal 1, just as the u values we have 
chosen in our case to avoid realizability problems. 
However, u is anisotropic in the TO case and u, is 


less than unity throughout the cloak region, 
tending to O at the inner boundary, which makes it 
difficult if not impossible to realize in practice. 
We will have further comments on this issue 
below, after we discuss the case of general 
transformation, which is non-linear in general. 


One can see from (20), which is valid when the 
transform from the free space to the transformed 
space is linear, that the relative permittivity in the 
latter space is proportional to the square of the 
ratio of the areas of the original to the transformed 
cells. If we follow this paradigm for the non-linear 
transform, we can obtain the parameter values in 
the transformed space by using a locally-linear 
approximation. As an example, we consider the 
geometry in Fig. (5), where the region of c < p < 
b is compressed into a region of c < p' < b via a 
non-linear transform, with the cell sizes varying 
exponentially in the radial direction from small to 
large in the transformed space, such that the cell 
size at its outer boundary 1s identical to that in free 
space. 

Let Ap be the cell size (uniform) in the 
original (free space) system in p-direction and let 
N, the number of rings in the original system, be 
given by 

N = (b—c)/Ap (24) 

In the transformed space, let the number of 
cells in the radial direction also be N, but let the 
cell size in this space increase exponentially from 
the inner to the outer boundary until it become 
equal to Ap (same as that in free space) at the 
outer boundary. (The advantage of this choice will 
be pointed out later). 

If we define T to be the multiplicative factor 
of a geometric series, we can express the cell sizes 
in the transformed domain as: 

Ap', = Ap: tt} (25) 
where i = 1~N , and i = 1 corresponds to the 
outermost layer. We can determine the value 
of T by using the relationship: 


N 
> _ 40", = 4p tht =b-a (26) 
i= 


which leads to the equation: 
1-1 b-a 
it ap (A1) 
Next, we take a cue from (20) to determine the 
relative permittivity €, by choosing: 
Eri =T OV (28) 





Let us consider a case example where the 
parameters of the geometry in Fig.5 are: c = 
1;b5=10;a=5 and t=0.8597 . Fig.5(c) 
shows the relative permittivity of the cloak 
derived by using (28). We can now see that the 
condition we imposed earlier on the cell size at the 
outermost layer of the transformed space, namely 
that it be equal to Ap there, leads to € |p'=p = 1. 


This in turn guarantees that the cloak region will 
transition into the free space without an 
impedance mismatch, which is obviously 
desirable. If we let c—>0, we can see that 
Eyy = TTY at the inner boundary will be large; 
in fact, we can show that it will tend to very high 
values as N > 00, 


At this point we return to Fig.I(c), where the 
results for the material parameters have been 
quoted from [5] for the same 2D problem, with 
identical polarization and geometrical parameters. 
It is interesting to note, first of all, that these 
medium parameters—both e€ and ym —are 
dramatically different from those given in Fig.1(b) 
for the same identical geometry, although those in 
Fig.1(b) are presumably valid for an arbitrary 
angle of incidence and polarization of the incident 
wave. This is a somewhat surprising, since one 
would expect from the laws of Physics that the 
medium parameters would go smoothly from the 
set in Fig.1(b) to those in Fig.1(c) as the incidence 
angle changes from oblique to normal, and not in 
a discontinuous manner as shown in these figures. 
The second observation we make is also quite 
interesting, namely that, as mentioned previously, 
the result for €, is identical to what we have found 
by using the WEB approach. Furthermore, lg 1s 


also equal fy in the TO case, which is again 
identical to our choice of u, which we have set to 
be isotropic at the outset. We point out that the 
material parameters in Fig.1(c) were obtained by 
applying the TO first, and then following it up 
with some simplifications after the fact (see [5] 
for details). 


Of course, it is evident that, despite the 
similarities, the WEB results are not the same as 
those from the TO, since the approximated TO 
results still call for an anisotropic u, for which 4, 
not only goes to 0 at r = a, but remains less than 
1 throughout the cloak region except at the outer 
boundary, where it approaches Ho. It goes without 


saying that such an anisotropic would be 
extremely difficult to realize in practice; 
furthermore, the cloak designed with the above 
parameters would fail to work when either the 
incidence angle is oblique or the polarization is 
arbitrary, even if we were successful in realizing 4 
values that are less than 1, and those that tend to 0. 
Hence such a cloak would be of limited practical 
use, not only because of the above shortcomings 
mentioned above, but also because it is likely to 
be very narrowband when realized by using 
Metamaterials. The fact that these observations 
are indeed true has been verified experimentally, 
and we mention this in support of the statement 
that we have given above. 


We next ask the question: “Is the WEB 
solution the complete answer to our cloak design 
problem?” After all, we have shown that the E, 
satisfies the same wave equations before and after 
the transformation from the p to p’ systems. A 
short answer is NO, because even though we have 
satisfied the wave equations in the two systems, 
we have not imposed the same boundary 
conditions at the outer boundary in the 
transformed system as those in free space, where 
the cloak transitions into free space with a 
discontinuity in the material parameters, which is 
bound to produce reflections from the interface, 
unless of course the outer boundary of the cloak 
recedes to infinity, which is unrealistic. As we 
well know, the solution to a differential equation, 
such as the wave equation, is not unique until the 
boundary conditions are stipulated, which 
explains why the WEB solution we have derived 
in the cloak domain is not the same as that in free 
space; after all, the two satisfy different boundary 
conditions at p = b. 


Though not shown here because of lack of 
space, we have also experimented with non-linear 
functions for mapping p into p’, in which the € in 
the p’ systems transitions from a relatively high 
value at the inner boundary, similar to the € 
profile shown in Fig.1(b), for the general TO 
cloak. For instance we have chosen e, = 275 at 
the innermost boundary, which transitions 
exponentially to €, = €ọ at the outer boundary to 
help mitigate the impedance matching problem at 
the interface with the free space, while we have 
still assumed that y = Up as we did for the case of 


linear transformation to circumvent the realization 
problem. Unfortunately, we find that this WEB 
design for the cloak also fails to deliver a 
satisfactory performance, leading us to conclude 
that neither the TO nor the WEB approach based 
on scaling provide workable “practical” solutions 
to the cloaking problem, for the reasons we have 
already given above that have to do with the 
realization problems, as well as with matching 
issues if we introduce approximations of the 
material parameters or sudden truncations that 
produce reflections at the outer boundary. With 
this background in mind, we will introduce a 
different performance metric for the cloak in 
Sec.3 in order to address this issue. 


Before closing this section, we mention that, 
for the general case where the mapping 
relationship takes the form: 


P' =f); P=6; z'=z (29) 


we define 

p = f) G0) 
where f~+ is the inverse function of f, and the 
Helmholtz equation of the E-field can be rewritten 
as 








dp 1 f- j E(P’, p) 
Op f- Eni Op M 
1 \*°07E,(p,¢) 
H) er 


+ ky“ E,(p',o) = 0 
(31) 


The above equation can be written in the 
form 








dp 1 “1p i E(P’, p) 
dp f- E Op 5 oe 
1 0°E, (0, $) 
. ol ab? 
+ W*Up€oEz(p',p) = 0 


(32) 


and there is no obvious way to choose the e€ in the 
cloak domain—as we were able to do in the case 
of linear transformation from p to p'—so that the 
cloaked target would perform as though it was 
invisible, provided the matching issue is resolved, 
of course. We have already mentioned that 
assuming an €, variation of the form which 
transitions from 275 at p = a, to 1 at p = b, also 
fails to provide the desired solution. 


Once the E-fields in the two domains have 
been matched, the corresponding H-fields which 
can be obtained from: 





I en 1 3Ez(p'r p) 
Hy, (Pp , $) — iW Up! ap (33) 
r p) = 2 980d) 
HCP P) g iW Ho dpr (34) 


will follow suit, and will match as well at the 
outer boundary if the cloak parameters tend to 
those of free space at the outer boundary, where it 
transitions into the external region without a 
mismatch. In fact, it is sufficient to match the E- 
fields only at the outer boundary, since we can 
invoke Huygen’s principle to argue that both E- 
and H-fields will match in the external region if 
we have done this. 


In the general case of the transformation 
from one co-ordinate system (€,, €,&3) to another, 
say (&1,&5,&3), the Laplacian takes the general 
form of 


VE = Xj Des ar Bo oe, (35) 
where 
Osa, Os, Ota 
0§, 06 0&3 
= sas “Cs: O82 
P= ety oe’ og, a8! = 
Osa Osa Oss 
0§, 06, 03 
and 
= 255 Oe, O55 OEE | 55 OE 
Pri = 38, og, | Be, O2 | OFs OF; eo 
for orthogonal systems, when k + j, 
Bg; = 0 (38) 
and the Helmholtz equation is written as, 
V'?E + k°E = “| : (JB =) 
J lð, X 0s 
J ð (e A 
0g, V e, 
+ (JB N+E = 0 
dg, V9 8, j 
(39) 


While there are no obvious ways we can use 
(39) in this case to determine the effective € 
values in the transformed (cloak) domain, we can 
make locally 2D approximations and apply the 
scaling relationship which we have derived above. 
Finally, we mention that if we are successful in 
rendering the two E -fields equal in the two 
domains at the outer boundary, then the 
corresponding H-fields in the external regions will 
follow suit, as we can show from Maxwell’s 
equations; hence, we argue that it is sufficient to 
work with either the E- or H-fields, whichever is 
convenient for the problem at hand. 


In summary, the WEB approach differs from the 
traditional TO approach in the following ways: (a) 
the WEB approach works with the Helmholtz 
form of the wave equation for either the Electric 
or the Magnetic field, rather than with Maxwell’s 
equations for both. Consequently, when we 
investigate incident waves of one polarization— 
exactly the same polarization investigated in a 
practical realization of the TO approach—the 
permittivity and permeability of the medium are 
isolated, and are required to have similar 
variations, typically anisotropic. In contrast to the 
TO equations, the permittivity and permeability 
appear as a product in the k term in the WEB 
approach, which enables us to choose non- 
magnetic materials and facilitates the practical 
realization of the cloak significantly; (b) the 
impedance matching property of the TO-dictated 
material parameters is rendered ineffective when 
the condition of practical realization is enforced 
and a reduced (approximate) form of the 
anisotropic material parameters is employed. 
However, in the WEB approach the boundary 
condition can be defined, without making any 
approximations, in a way such that the 
permittivity of the cloak medium is matched at the 
cloak-medium-air interface, in order to ensure that 
there is no reflection at the boundary. 


III. Conformal Mapping Approach 


In this section we will discuss an alternate 
form of WEB approach, which is based on the 
Schwartz-Christoffel type of conformal mapping, 
and which provides a systematic way, not only to 
transform the wave equation, but also to impose 


the condition that the permittivity approaches €p at 
the outer boundary. 


It should be pointed out that if we choose 
Ho = 1, as we have done in the present case, 
the magnetic fields will not satisfy the mapping 
relationship, even if the electric fields do. Hence 
we would need to ensure that the H-field in the 
transformed domain will tend to the H-field in the 
original space as the observation point moves to 
the outer boundary of the truncated cloak. To 
satisfy this requirement, we should not only focus 
on the transformation of the wave equation, but 
should also pay a close attention to the boundary 
condition. Below we show how we can use the 
Schwartz-Christoffel (SC) type of conformal 
mapping to ensure the impedance matching as 
well as the continuity of the fields at the boundary 
by appropriately choosing some of the transform 
parameters in the SC mapping. 


Although the conformal mapping approach to 
transforming the wave equation is strictly limited 
to two-dimensional (2D) problems, we can again 
gain much insight from this example and draw 
some conclusions that are also applicable to the 
three-dimensional case. As we pointed out earlier, 
similar two-dimensional cases, e.g., bumps in 
PEC planes have been researched earlier [29] by 
using the TO algorithm, and carpet cloaks are but 
one of these examples. 


We consider a triangular bump in a PEC plane, 
as shown in Fig.6(b), which is defined by the 
coordinates (u,v), so that we can transform it 
conveniently to or form a flat plane in the (x,y) 
system, as in Fig.6(a), by using the well-known 
Schwartz-Christoffel transformation [30, 31]. 


We assume, in common with the existing 
literature, that the fields are independent of the 
longitudinal direction z, and that all the fields can 
be derived from the E,-component of the Electric 
field, which is the only non-zero component of the 
Electric field in (4). 


In common with Sec.2 we consider the TM, 
case once again, and set E = ZE,. Next, as before, 
we eliminate H from the curl equations in (4) to 
arrive at (9), where we have assumed that u is 
simply fp in anticipation of the fact that we would 
search for dielectric-only materials when we carry 
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Fig.6. (a) Free space; (b) physical space 


out the coordinate transformation of the 
Helmholtz equation (wave equation). 


Our next step is to conformally map the 
(x,y) domain into the (u, v) system by using the 
Schwartz-Christoffel (SC) transformation, such 
that the flat plane in the ¢ = (x + jy) co ordinate 
system (free-space domain) maps into the 
triangular bump in the w = (u + jv) coordinate 
system (physical space and cloak domain), and 
vice versa. The SC transform relationship reads: 


Ç 
ose | (E + k)-1/3(€)2/3( — k)-1/3dE + C, 
0 


(40) 


where the parameters C,k,C,, are yet to be 
determined. A systematic approach for doing this 
will now be presented below. 


From Fig.6 we see that the three points 
€, =—k,& =0,ţ%;}=k in free space are 
transformed to the corresponding points w, = 
—1,w, = V3j,w3 =1 in the physical space. 
Letting (, = 0 and œ, = V3j in (40), we can 
readily obtain C, = V3j. In view of the symmetry, 
either of the points 1 and 3 can be chosen to 
determine the values of k and C. Using ¢3 = 
and letting w3 = 1 in (40), we get: 


ciel f on r| df = 267 sz (41) 


By taking the modulus of each side of (41), 
we obtain: 


K &2/3 
ci f helaa (42) 


The choice of k in (42) appears to be 
arbitrary at this point if we do not impose any 


condition on |C|. And, in fact, if we use the SC- 
Tool of MATLAB code [32], it makes the default 
choice k = 1, which leads to |C| = 2.3192. 


However, we will soon see that this choice of 
k is not suitable for our purpose, since the 
asymptotic behavior of |dw/d¢| is not what we 
desire for it to be. To explain why this is the case 
we insert a box-shaped boundary in the (u, v)- 
domain as shown by dotted line in Fig.6. Since, 
from a practical point of view we must truncate 
the cloak with which we are going to cover the 
triangular bump, we define the outer boundary of 
the cloak so that it coincides with the surface of 
the box. We will show below that by introducing 
appropriate dielectric-only materials in the cloak 
region of the (u, v)-plane, we can render the E,- 
field in the transformed (cloak) domain to be the 
same as that in the half-space of the (x, y)-plane, 
which is free space. To show this, and to 
determine the material parameters of the cloak 
region, we turn to the Helmholtz equations of (9) 
in the free space of ¢ domain rewrite as: 


Vey? (x,y) + kõ¥ (x,y) = 0 (43) 


In the w-domain, the medium parameter are 
Ho and €;,€ 9, where €, is yet to be determined. 
Then ¥(u, v) satisfies: 


V2, (u,v) + ki Y (u,v) = 0 (44) 
where 
ki = kon? = (2n f)’ uoEoEr = koe, (45) 
We can transform (44) into: 
Vey (x,y) + kĝer(d w/d ()*¥ (x,y) = 0 
(46) 
To render (46) equal to (43), we let 


Er = (dw/d¢)~* (47) 


This is the choice for the material parameters 
we need to make in the w (cloak-domain), so that 
the field solution in this domain is identical to the 
one in the ¢ -domain, ie, P(x,y), which 
comprises only the specularly ”reflected” field 
from the PEC plane, and which implies that the 
scattered field is zero. We can then argue that with 
the choice of (47) for €,, the cloaked bump on a 
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Fig.7. The mesh in the original (left) and the physical 
space/transformed space (right). 
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Fig.8. The relative permittivity of cloak. The color bar 
is not shown since €, becomes infinite at the tip of the 
bump. 


PEC plane would also produce no scattered field, 
but only the specularly reflected field; hence the 
bump would be invisible. 


Since | dw/d¢ |o» = |C|, we must choose 
|C| = 1 if the boundary conditions at infinity are 
to be the same. This implies, in turn, that 
Erlu»» = 1. Furthermore, if the region external 
to the “carpet” is to be free space so that the cloak 
can be truncated, the permittivity of the carpet 
should approach the permittivity of free space to 
minimize any reflection at the interface. This also 
requires that we enforce the condition |C| = 1. 


Now we can understand why we mentioned 
earlier that the choice of |C| = 2.3192 is not good 
for designing a cloak with minimum reflection at 
the interface without the need of a matching layer, 
which we required for the case example presented 
earlier in Sec.2 for cloaking the cylindrical object 
designed by using the WEB approach with a 
linear transformation relationship between p and 


p . Fortunately, for the case of conformal mapping 
we can adjust the value of |C| by choosing an 
appropriate k. For instance, using (40), we can 
show that if k = 2.3192, |C|=1, and the 
truncation condition €,|,,,. =1 is indeed 
satisfied. 


Figs.7(a) and (b) display the meshes in original 
and the transform spaces. Fig.8 presents the 
material (€ variation) that we need to realize for 
the carpet. We note that €, — œ at the tip of the 
triangular bump, and is equal to zero at the two 
points u = +1,v = 0, where the triangular bump 
joins up with the PEC plane. It is obvious that 
these parameter values would be difficult to 
realize, even though unlike TO the required e€, 
values are isotropic in this case. Even if we were 
to use a bump with a smooth curved profile to 
replace the triangular bump, it will still require 
€, <1, which is difficult to realize. 


Obviously, strictly speaking, such a cloak is 
not physically realizable, not only at the tip region, 
but where the bump joins up with the PEC plane 
and the required € — 0. 


In any case, we have seen that even for the 
example of the simple cylindrical geometry 
presented above, the material parameters required 
for the cloak designs based on TO are difficult to 
realize, whether they are exact (see Fig. 1(b)) or 
approximate (Fig.1(c)). Additionally, we note that 
at p’ = b, which is the outer boundary of the cloak, 
there can be a discontinuity in the material 
parameters, if we make approximations such as 
assume U = Lg, aS we did in the WEB design, or 
if we choose simplified parameters for € in TO or 
WEB designs for the sake of realization. This is 
because they might produce untoward reflections 
of the incident wave impinging upon the cloaked 
cylinder from outside of the cloaked region. 


In summary, there are two major problems 
with both the TO and WEB approaches to cloak 
design, where we seek to reduce the radar 
signature of the target in all directions (bistatic), 
and to make it invisible in the limiting (ideal) case. 
First of these is the material realization problem, 
despite the fact that we only need isotropic and 
homogeneous dielectric materials in the present 
WEB-design, as opposed to inhomogeneous and 
highly anisotropic materials called for in the TO- 


design, as we have pointed out before. The second 
problem is associated with the boundary; whose 
choice appears to be arbitrary. It is claimed in TO, 
that this problem is circumvented by requiring 
€ = Lt; however, that solution is hardly realistic 
from a practical point of view, since realizing 
even isotropic values with the required range of 
parameters is a formidable if not impossible task 
for TO-based designs, and the situation with the 
realization of € is not much better either. This is 
evidenced by the fact that various approximations 
have been suggested [5] to get around this 
difficult issue—after the initial design has been 
found to be fraught with realization problems— 
albeit at the expense of compromising the 
performance of the cloak. In principle, we could 
address the matching issue in the WEB design, but 
the realization issue would still persist. For this 
reason, in Sec.3 we will propose a different 
performance metric for the cloak, which is widely 
used by the microwave community, when we 
discuss the Scattering matrix approach to 
formulating the problem at hand. 


We will close this section by providing one more 
example which underscores the importance of 
satisfying the boundary condition at the interface 
between the transformed and external regions. 
Towards this end, let us consider the example 
shown in Fig.9 to illustrate our point. On the left 
side of this figure we see a sectoral horn whose 
aperture size is a, Let us say we map the geometry 
of this horn to the one on the right whose aperture 
dimension is shorter, say b, where b < a. Then, if 
we were to conjecture that we can just transform 
the inner regions (shaded), and find materials that 
make the fields in the two shaded domains the 
same, say by using either the TO or WEB 
algorithms, we would think that we could get the 
horn with the small aperture (right) to perform just 
the same as the one with a larger aperture (left) 
and to provide the same gain values. While it 
would be very desirable to reduce the aperture 
size of an antenna and achieve super gain 
characteristics without compromising its 
performance by using the TO algorithm, 
unfortunately this turns out not to be possible. 
This is because to ensure the field equivalence 
between the two domains we must transform the 
entire region in which the horn is embedded, and 
not just the shaded region. The material 
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Fig.9. Sectoral horn in the (a) original system with 
aperture size a and (b) transformed system with a 
smaller aperture size b. 


parameters of the horn on the right would extend 
well into its external region and will not be 
confined to the interior (shaded) region alone, in 
order to correctly map the fields from the left to 
the right figures, or vice versa. Consequently, the 
true aperture size of the horn in the right will in 
general be larger than the opening of its mouth, 
and will depend on where we place the truncation 
boundary of the loading materials. Note that, once 
again, the boundary condition will play an 
important role, and we will face the problem of 
reflection from the interface if the medium 
parameters do not tend to the free space values, 
namely €o an Up, at the outer boundary, or the 
impedance is matched some other way, e.g., by 
using anisotropic materials with parameters that 
are difficult to realize in practice. 


IV. Generalized Scattering Matrix 
Approach 


In the previous sections we have discussed 
two different WEB  (wave-equation-based) 
approaches for designing cloaks and carpets that 
can reduce the scattering from certain types of 
objects, or even render them invisible, by using 
dielectric materials alone. We have found that 
material (relative permittivity) parameters 
required to accomplish this can vary over a wide 
range; furthermore, in the TO designs, €, and u» 
can either be much smaller than 1, much greater 
than 1, or combinations thereof, depending on the 
geometry of the object . We point out that this 
problem is compounded in the TO-based 
approaches because there both € and u are tensors 
in general, and this makes it even harder to realize 
them, even using Metamaterials. It is not 
uncommon therefore, to find designers 


introducing approximations and replacing the TO- 
dictated materials with much simpler ones (see [5] 
for instance) that deviate substantially from the 
TO-specified ones in (48), so that they are more 
realistic in terms of realization potential than are 
the original TO ones expressed in (49), albeit at 
the expense of compromising the performance in 
comparison to that achievable by using the ideal 
set of parameters. In [5] these modified 
parameters are: 
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which are significantly different from the original 
ones dictated by the TO, and are plotted in 
Fig.1(b). 


The root cause of this difficulty can be linked 
to the fact that in the TO-paradigm we insist that 
the field quantities in the original and transformed 
systems be identical, as we do when we are 
locked into the TO-based paradigm. In this 
Section we propose a different performance 
metric in the context of Generalized Scattering 
Matrix (GSM), which we introduce below. We 
show that the TO prescription of field equivalence 
between the original and transformed systems is a 
special case and that the GSM provides us more 
flexibility than does the TO, by modifying the 
performance metrics of the electromagnetic device 
that we are designing, dramatically increasing its 
chances for practical realization by using available 
materials. Thus although the performance of the 
device following the GSM prescription will be 
less than ideal, it may be quite satisfactory from 
the point of view of practical applications that we 
may have in mind. We will provide several 
examples to illustrate this by considering RCS- 
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Fig.10. Generalized Scattering Matrix approach in the 
context of the Field Transformation (FT) method 


reducing blankets, flat lenses, reflect arrays, etc., 
all of which face the realization problems when 
using the TO approach, which calls for anisotropic 
materials (both € and u), and whose realization 
has been elusive at best. 


To introduce the GSM approach [26, 33] in the 
context of the so-called Field Transformation 
method [34], we refer the reader to Fig.10, where 
we have defined the input and output ports to 
correspond to interfaces that bound an 
electromagnetic device. The field distribution in 
the input port, which is illuminated by the source 
located at the left of the port, can be expressed in 
terms of a set of coefficients al (vector) 
associated with the basis functions employed to 
represent this “incident” field in the absence of the 
device when there are no reflections. A 
convenient choice for the set of basis functions is 
the discrete plane wave spectrum when we are 
dealing with objects in free space, where the 
discrete plane waves play a role similar to those of 
the set of waveguide modes in “closed” guided 
wave structures. 


Next, we insert the electromagnetic device, 
whose Scattering Matrix we desire to describe, 
inside the region bracketed by the input and 
output ports. We define a set of coefficients bt, 
again associated with the same basis functions as 
we used to define a}, to represent the outgoing 
fields scattered by the electromagnetic device, i.e., 
the “reflected” fields that originate from the 
device and propagate back towards the source. We 
can similarly define a set of coefficients 
c} associated with the field distribution in the 
output port, through which these fields propagate 
in the free-space region to the right of this port, 
and are termed the “transmitted” fields. Our next 


step is to place the illuminating source to the right 
of the output port, which we have previously 
defined when the source was at the left, and 
reverse the roles of the input and output ports to 
correspond to the new source location. The 
incident, reflected, and transmitted fields are now 
characterized by a new set of coefficients a?, b?, 
and c?, where c? fields now propagate to the left 
of the device, whereas the b? fields do the 
opposite, i.e., propagate to the right. 


Below we present several examples of 
practical devices where we use the GSM to define 
their performance metric. 


We consider a radar target whose RCS we 
wish to reduce by covering it with a suitable 
material, which we can either fabricate in the Lab 
or which is commercially available, so we do not 
have to wrestle with the feasibility issue in the 
design process later. The performance metric we 
define is that the “magnitudes” of the elements of 
S11, comprising of the scattered fields from the 
target in the left hemisphere be below a certain 
threshold level. It is important to note that we do 
not impose any condition on the fields scattered in 
the forward direction, i.e., on S42, either its 
magnitude or its phase, neither do we stipulate any 
condition on the phase behavior of the elements of 
S141. This is in contrast to TO, which aims for the 
ideal situation where the magnitudes and the 
phases of both $3, and S12 at the input and output 
ports, respectively, be identical to the same as 
they would be if the target were invisible, or was 
absent altogether, implying that the fields 
scattered by the target to the left of the input port 
or to the output port be all identically zero. 


We are now ready to define the scattering 
matrix [|S], via (50) below, which we will use to 
characterize the device, as follows: 


b = [S]a (50) 
or explicitly, 
S11 ra 
S| = 51 
[S] Sa S» (51) 


where b = {bt,b?} and a = {at, a°} represent 
the weights of the outgoing and incoming field 
representations at the input and output ports, 
respectively. Equation (51) provides us a 


convenient way to characterize an electromagnetic 
device in terms of its response to a plane wave, 
regardless of whether the illuminating source is 
incident from the left or the right of the device. 


It is not difficult to see that the GSM metric 
is considerably easier to realize than the 
corresponding TO-dictated metric we just 
mentioned. Thus if in the radar scenario it is 
deemed sufficient to be concerned only with the 
backscattering region, the GSM metric will suffice, 
and we can use available [35-37] electric and 
magnetic materials to design a thin coating for the 
target (typically only a few millimeters thick) 
which will reduce the scattering (RCS) over a 
very wide frequency band, and for arbitrary angle 
of incidence and polarization of the radar signal 
impinging on the target, which is highly desirable 
in practical radar applications. The multilayer 
coating essentially serves as a matching region 
which transforms the surface impedance of the 
target (zero for PEC scatterers) to the wave 
impedance of free space at the interface between 
the coting and external region, which is free space 
by choice because we wish to truncate the cloak, 
which is theoretically infinite. The transition 
between the coating and the external region is thus 
handled naturally in this scheme in one-step 
process. 


We should mention here that attempts have 
also been made to make the coating thin for TO- 
based designs, albeit for simple cylindrical-shaped 
targets of infinite length, for the case of radar 
signal at normal incidence (0 = 90°). We present 
the results of numerical simulations below that 
show that the performance of three-layer 
discretization of the thin cloak is not very 
satisfactory. 


We examine the problem of reducing the 
level of scattering from radar targets—a problem 
that has been extensively researched into in the 
context of the TO. Once again, we will cast the 
problem in the language of Scattering Matrices, to 
help us understand why the TO approach leads us 
to untenable situations and/or to solutions which 
call for MTMs that have many drawbacks as we 
have pointed out earlier. We will also show how 
the FT approach mitigates the problems alluded to 
above by restating the design objectives and 
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Fig.11. Scattering from an arbitrarily shaped target 
described in the context of Generalized Scattering 
Matrix (GSM) method in the physical domain (a) and 
virtual domain (b). 


modifying them slightly from those associated 
with the TO-based designs. 


Let us consider an arbitrarily shaped radar 
target placed in the region between the input and 
output ports, as shown in Fig.11. Next, let us 
suppose that our objective is to reduce the level of 
scattering from the target, both in the forward and 
backward directions to the extent that the target 
becomes totally invisible to the incident field, say 
from an interrogating radar. We could cast this 
objective in the language of the Scattering 
matrices, by specifying that $,, be identically 
zero at the input port and S42 be such that the 
field distribution at the output port is identical to 
the incident field, as though the scatterer was 
totally transparent or invisible. 


To circumvent the problems that we 
encounter when attempting to use the TO in order 
to meet the ideal but unrealistic goal of making 
the target altogether invisible, we turn to the FT 
(Field Transformation) approach and modify our 
stated objectives, again in the context of S- 
parameters. In contrast to the TO, this time we ask 
that $11 be small—in terms of magnitude only— 
but not 0, as we demanded in the case of TO. 
Furthermore, we do not impose any restrictions on 
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Fig.12. Scattering from an FT-treated arbitrarily shaped 
target described in the context of Generalized 
Scattering Matrix (GSM) method. 


S12, aS we did in the case of the TO when we 
stipulated that the scattered fields at the output 
port be identically zero, so that the total field there 
be just the incident field in the absence of the 
target. While we concede that the performance of 
the FT-based cloak or the blanket won’t be as 
ideal as the performance of the TO-based design 
would have been if we could realize it, we 
certainly stand to gain considerably when we 
follow this strategy, since we can now obtain 
realizable solutions for the cloak that are very 
wideband and can cover the entire frequency 
range of 2-18 GHz, for instance, if we so desire. 
Furthermore, the cloak (or blanket) can now be 
very thin (only a few millimeters) and it would 
work for arbitrary incident angles and 
polarizations as well. In short, we can say that this 
strategy of following the FT-based algorithm to 
design the cloak enables us to circumvent all the 
problems we encounter when employing the TO- 
based design strategy instead. Of course, we 
compromise the performance of the cloak in the 
forward scattering direction when we use the FT- 
based strategy, though that is not a problem for 
either monostatic or bistatic radars that are only 
concerned with the backscattered fields. 


Some examples of such realistic blanket 
designs based on the FT strategy can be found in 
[38] and are omitted here for the sake of brevity. 
The above work also shows how to design the 
blankets for arbitrary targets by borrowing some 
ideas from the TO algorithm to improve their 
design. 


Before closing this discussion, we would like 
to mention that the strategy for designing carpet 
cloaks is very similar. We can use the image 
theory to recast the problem of designing a carpet 
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Fig.13. Real and imaginary parts of the two types of 
absorbing materials used in the multi-layer absorbing 
blanket design: (a) permittivity and (b) permeability. 


cloak so that it totally equivalent to the cloak 
design problem for radar targets that we have 
discussed above in this section. 


As is well known, the RAMs (radar 
absorbing materials) have been around for a very 
long time, some for many decades, dating back to 
when stealth aircrafts came into vogue in the 
sixties, and earliest theoretical and experimental 
work starting around 1930s. We realize that 
information on some of these RAM materials is 
not available openly because of their “classified” 
or “secret” nature, understandably so because they 
are used in military applications to design stealth 
aircrafts and missiles. Nonetheless a plethora of 
information about similar absorbing materials is 
available in the open literature, including the 
details of their fabrication, which have been 
described in [35, 36], for instance. 


Here we will use two different types of 
materials namely CoFe Nano-Flakes (NF) and 
CoFe Nano-Particles (NP), whose frequency 
variations are shown in Fig.13. We point out that 
these materials can be realized with relative ease, 
as is evident from [35, 36], where the details of 
their fabrication can be found. 
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Fig.14. Reflection coefficient of multilayer absorber 
backed by a PEC plate for different number of layers. 


We note that the € and u values of these 
materials are both complex and dispersive. We 
hasten to point out, however, that unlike the 
MTMs, which typically have variations with 
frequency because of resonant inclusions they use 
to achieve effective € and u values that are very 
high, very low, or even negative, the material 
parameters of the NP and NF materials vary 
relatively smoothly with frequency, and this is 
crucial for realizing wideband performance, as we 
will soon see. 


To illustrate the fact that we can indeed 
achieve wideband performance in terms of 
reflection reduction over a wide frequency band 
with relative small thicknesses of 2, 4, 6 and 7 
layer absorbers we refer to Fig.14, in which a 
10dB (or better) reduction in the reflection 
coefficient is presented. Although not shown here, 
the results for the reflection coefficient reduction 
are also satisfactory when either the polarization, 
or the incident angle is varied and this is also true 
when both are changed simultaneously. 


We now move to the second step in our 
design procedure, which is to adapt the blanket 
designed for the infinite PEC plane to an 
arbitrarily shaped object. Initially we consider an 
object with a smooth surface whose radius of 
curvature is moderate-to-large everywhere. We 
will generalize the procedure in the third step, 
using the principles of the TO when the above 
assumption regarding the smoothness of the object 
is not valid, as for instance when the object has 
sharp edges or bumps, as a general target would in 
practice. 


When the object has a relatively smooth 
geometry, we initially wrap the multilayer 
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Fig.15. Back-scattering RCS of PEC objects covered 
with 2 and 7-layer absorbers: (a) plate; (b) pyramid and 
(c) cylinder. 


absorbing blanket, which we have designed earlier 
for the planar surface around the PEC target 
whose scattering cross-section we are attempting 
to reduce, and test the effectiveness of the blanket 
for the new object. For a wide variety of targets 
that we have examined, a number of which are 
shown in Fig.15, we have found that the blanket 
does reduce the monostatic as well as the bistatic 
radar cross-section in the “reflection” region near 
the surface of the object for different angles of 
incidence and polarizations of the incoming wave. 
The results for a two-layer absorber are shown in 
Figs.16 and 17 for a rectangular cylinder of finite 
length, which we have studied as a test case. 


A simple test, which is typically applied to 
cloak designs, is to examine the wavefront of the 
total (ancident+scattered) field, and see how the 
level of distortion of the wavefront decreases 
when the scatterer is covered by the layered 
absorbing coating. We present the plots of these 
wavefronts of the total fields for normal and 
oblique incidence cases for both polarizations in 
Figs.16 and 17, respectively. 


We observe that the object, which is a finite 
cylinder of height 18cm, generates distorted 
wavefronts owing to the contribution of the 
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Fig.16. Phase behavior of the E-field near the 
rectangular PEC cylinder, which is wrapped around by 
an absorber blanket, normally incident on the cylinder. 
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Fig.17. Phase behaviors of the E-field near the 
rectangular PEC cylinder, which is wrapped around by 
an absorber blanket, for an obliquely incident plane 
wave. 


scattered field from the object, even when covered 
by a two-layer blanket, designed for the infinite 
planar PEC object, for the nominal frequency 
range of 4.6-18 GHz, with a nominal reflection 
coefficient of -10dB or less. However, we also 
note from Fig.16 that the distortion in the phase 
front is relatively small once we go above the low- 
end of the design frequency range, viz., 4.6 GHz 
for the planar geometry, confirming that the 
planar design performs reasonably well even 
though we are dealing with a rectangular cylinder 
now. We hasten to point out that the results 
presented in Fig.16 are not for a coating which has 
been optimized for the object at hand, and we 
expect some compromise in the performance of 
the coating. However, we can improve this 
performance by optimizing the parameters of the 





Fig.18. Target detection in conventional radar scenarios: 


mono-static and bi-static schemes. 


two-layer design, specifically the relative 
thicknesses of the two-layer, even as we maintain 
the total thickness intact. We expect the changes 
to be relatively minor, however, except for the 
corner regions and, hence, the optimization 
process should be realistic as well as numerically 
feasible. The above remarks are also applicable to 
the oblique incidence case, for which some 
sample results are presented in Fig.17. 


It is important to point out that the strategy 
for designing the absorptive coating, presented 
herein, is very different from that employed for 
ideal traditional TO cloak, since the latter is 
designed to render the (object+cloak) composite 
to have a zero scattering cross-section in all 
directions, whereas the blanket design introduced 
here seeks to reduce both backscattering and bi- 
static scattering scenarios but only in the 
reflection region, and not in the forward-scattering 
direction. We hasten to point out, however, that 
this type of performance is perfectly well suited 
for modern radar systems, symbolically depicted 
in Fig.18 above, where only the scattering in the 
reflection region is of concern. 


For the final step, we consider the problem of 
absorber design when a shape perturbation is 
introduced in an object. Let us say that our 
modified target 1s the same rectangular cylinder 
we just considered above, except for a bump on 
the top surface. The extra corners introduced by 
the perturbation, be they smooth or sharp, would 
obviously introduce additional distortions in the 
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Fig.19. (a) Original 2-layer absorber wrapped around a 
rectangular cylinder with shape perturbation (left); (b) 
2-layer absorber with TO-modified material properties 
around regions with the shape perturbation (left); (c) 2- 
layer absorber with TO modified material properties 
for all regions (left) and 2-layer absorber wrapped 
around a rectangular cylinder (right). 


planar phase front, and potentially increase the 
scattering level. Our objective here is to restore 
the field behavior so that it is close to that of the 
original object that we had prior to the 
introduction of the perturbation. 


We now outline the procedure for the blanket 
design for the new object, shown on the left in 
Fig.19, i.e., Fig.19(1), which is in the physical 
domain, and is a modified version of the one 
shown in the right side of the same figure; i.e., 
Fig.19(2), which corresponds to the virtual 
domain. Note that unlike the circular cylinder 
example we discussed earlier, the medium in the 
virtual domain, surrounding the object, is no 
longer free-space, as was the case shown in Fig.18. 
Note also that the dimensions of the objects in the 
two domains are comparable, and are totally 
different from the legacy TO-design case, in 
which the scale factor between the dimensions of 
the object in the physical and virtual domains 
tends to infinity to render the target invisible. 


The field behaviors for the perturbed object 
with locally modified material parameters can be 
seen from Figs.20 and 21. Fig.20 shows that the 
amplitude of the scattered E-field 1s reduced, and 
that the phase front of the total E-field is 
approximately restored as well in the reflection 
region. Fig.21 compares the amplitudes of the 
scattered electric fields for five different scenarios 
listed in the figure. 
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Fig.20. Phase behavior of scattered E-field for (a) 
perturbed object wrapped by blanket with original 
medium parameter; (b) perturbed object wrapped by 
blanket with locally modified medium parameter. 
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Fig.21. Comparison of the amplitudes of electric fields 
scattered by different objects. 


We should clarify the fact that although we 
are referring to this geometry as a slab, what we 
are really dealing with is a wide rectangular 
cylinder, with a small thickness. 


We note that the modified slab does 
introduce additional scattering, and that the 
absorber does help reduce the same. We also note 
that the modified absorber improves the 
performance over the initial one, but only slightly, 
which shows that the planar version of the cloak is 
not all that inferior to the one modified for this 
type of geometry. Additional optimization of the 
modified cloak is expected to improve the 
performance even further, 1f so desired. 


We look at an important application of RFI 
reduction, i.e., to mitigate the problem of antenna 
blockage in a shared-platform environment, as 
shown in Fig.22, in which the introduction of the 
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Fig.22. Model schematics for the blockage problem: 
(left) side view and (right) isometric view of the victim 
antenna (parabolic dish) and the aggressor antenna 
(monopole) sharing the same platform. 
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Fig.23. Radiation patterns for the antenna blockage 
problem: (blue) dish antenna, (green) dish antenna and 
the monopole antenna in the vicinity and (red) dish 
antenna with absorber-treated monopole antenna. 


aggressor monopole antenna can raise the far-end 
sidelobe levels of the parabolic dish significantly. 
To mitigate this effect, we can wrap the monopole 
by using a multilayer absorber. Alternatively, we 
can use a discretized cloak structure to cover the 
monopole. Fig.23 shows the far-field patterns of a 
parabolic dish antenna—the victim antenna— 
when the aggressor antenna, is placed in the 
vicinity of the dish. The patterns of the 
dish/monopole composite are also included in the 
figure for comparison with the two treatment 
plans applied to the aggressor antenna. We note 
that the introduction of the absorber treatment 
reduces the sidelobe levels of the dish antenna, as 
compared to the case for the far-end antenna 
combination without the treatment. We also note 


that an elaborate cloak design, however, does not 
mitigate the pattern interference problem, or rather, 
aggravates it. 


V. Metasurface and Lens designs 


The subject of field manipulation is a topic of 
great interest today, as evidenced by a rather large 
number of publications on this topic in recent 
years. The present authors have discussed this 
topic elsewhere in a recent publication [33] in the 
context of the Generalized Scattering Matrix 
approach, which we have described above in 
Sec.4. Here we simply summarize the approach 
briefly and mention some representative examples 
for the sake of completeness. 


We begin with the problem of designing a 
flat surface that mimics a parabolic reflector, 
which may be either symmetric or offset-fed. We 
assume that the reflector is illuminated in the 
usual manner, by using a feed horn which 
emanates a spherical wave from its phase center 
(see Fig.24). The reflector serves to convert this 
wave front to a planar one to generate a directive 
beam in a specified direction, which depends on 
the choice of the feed location. 


The TO approach to handling the problem of 
designing the flat reflector is relatively 
straightforward: Transform the parabolic reflector 
geometry into a planar one and embed it in a 
medium that surrounds it so as that the fields 
reflected off the surface of the planar reflector 
covered by the TO-designed coating mimic those 
associated with the original parabolic surface. The 
Jacobians of the transformation provide a 
relatively easy way to determine the parameters of 
the medium in which we must embed the flat 
reflector so that it delivers a performance that is 
identical to that of the parabolic reflector. But we 
encounter similar types of realization problems in 
the TO design of these surfaces, as we did when 
realizing cloaks, namely finding the desired low- 
loss, low-dispersive and wideband materials with 
which to cover the surface of the flat reflector. We 
note, however, that such surfaces that are also 
known as reflectarrays, are almost always 
designed to match only the phase of the reflected 
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Fig.24. Radiation from two antennas described in the 
context of Generalized Scattering Matrix (GSM) 
method: (a) parabolic dish antenna and (b) reflectarray. 


field, rather than to control the amplitude 
distribution of this field in its equivalent aperture 
as well. Thus, in the context of S-parameters, we 
can simply use a one-port S-parameter description 
(see Fig.6) since the surface is backed by a PEC 
plane, and require only the phase of $44, i.e., that 
of the reflected field, to satisfy the desired 
specifications. We can then design the surface 
with isotropic and dielectric-only materials by 
locally controlling just the phase of the reflection 
coefficient. In this approach we can circumvent 
the need to use both € and u, which would not 
only be anisotropic but would require the use of 
Metamaterials as well, if we follow the TO 
paradigm strictly and proceed to match the fields 
in the original and transform domains in their 
entirety. 


Similarly, for the case of a flat lens which 
mimics a conventional convex lens, we impose 
the condition only on the phase of S42, where we 
locate the input and output aperture planes just to 
the left and right of the aperture, respectively (see 
Fig.25). Here again, we can wuse isotropic 
dielectric-only materials and avoid the problems 
typically with the TO designs. We do not impose 
any condition on either the magnitude of $44, or 
on the magnitude of S32, though we would 
obviously wish to minimize S41 for an efficient 
design. 


For the last example we consider the problem 
of manipulating the amplitude distribution in the 
aperture of a feed horn to control its edge taper. 
While the primary use of the TO algorithm is to 
relate the material parameters of two transformed 
domains to each other such that the field solutions 
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Fig.25. Lens problems described in the context of 
Generalized Scattering Matrix (GSM) method. 


in these two domains are identical, it is possible to 
generalize it to determine the material parameters 
of the intervening medium that would transform 
the field from the input to the output port (see 
Fig.10 in Sec.4) in a desired manner. As an 
example, we may wish to modify the aperture 
distribution of a horn antenna (see Fig.26) to 
achieve certain characteristics of its radiation 
pattern, e.g., the sidelobe level, which is 
determined by its edge taper of the distribution. 
We will now describe how we might control this 
edge taper by inserting a material slab in the 
aperture of the horn, or just above it. We will now 


discuss two different ways to accomplish this goal, 


namely to use the TO and Ray Optics approaches. 


In the past, the TO has been used to control a 
beam in several ways, e.g., to bend and expand the 
beam [39], to shift and split [40], and to generate 
multiple beams [41]. In this section we will 
examine the question whether or not we can 
transform a given field distribution at the input 
port to the desired one at the output port by using 
the TO. As shown in [42], the field can be 
concentrated to have a higher energy density in a 
region by using appropriate material parameters 
dictated by the TO. In contrast to [42], we attempt 
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Fig.26. Magnitudes of aperture distributions: original 
(blue) and desired (red). 
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Fig.27. Horn aperture with a material slab to control 
the aperture distribution. 


to use a method that manipulates the field 
distribution profile by letting the fields propagate 
through a region in which a uniform mesh is 
modified to a denser or coarser mesh--depending 
upon whether we want a higher or lower 
amplitude in the output plane--and the cells are 
filled with media whose parameters are dictated 
by the TO. It should be pointed out that the mesh 
inside the slab (see Fig.27) should be carefully 
designed to ensure a smooth transition from the 
input to the output aperture. 


AS mentioned earlier, the field 
transformation via TO would typically call for the 
materials in the slab to be anisotropic. Even if we 
could realize the anisotropic lens designed by 
using the TO, we should recognize the fact that 
the performance of the lens would not be as we 
would desire for it to be. This is because we 
would have transformed only the fields inside the 
lens, and the fields outside the lens have not be 


transformed. As a result, there would be 
reflections at the boundaries of the lens because 
we have ignore the boundary conditions at the 
interfaces. 


Before closing it is worthwhile mentioning, 
that the lenses and reflectarrays designed by using 
the present S-parameter-based approach are not 
only easier to realize in practice, but they have 
been found to exhibit superior performances as 
well when compared to TO-based designs. 


It is worthwhile to point out it is possible to also 
achieve the desired field manipulation by utilizing 
isotropic dielectric materials only, which is highly 

desirable from realization point of view. The 
design of the slab, which acts somewhat like a flat 
concave lens in this case, has been carried out in 

[29] by using a ray-optical approach in the same 

manner as that employed for the flat lens design 

we have described above. Fig.28 shows an 
example of this type of design, where the rings 

with varying permittivity are used to construct a 
flat concave lens that covers the horn. The relative 
permittivity values of the rings, listed in Table 1, 

is seen to increase form center to the edge. This, 
in turn, introduces a larger phase delay at the edge 
than at the center, which results in an expansion of 

beam as shown in Fig.28(c), where the aperture 
distributions with and without the lens are 
presented. Note that, in contrast to the TO-based 
design, the medium in the present design is 
dielectric-only, is isotropic and requires only 
realizable materials rather than Metamaterials. 


VI. Conclusions 


In this work we have investigated two 
different Wave-equation-based (WEB) approaches, 
namely coordinate transformation and Conformal 
mapping, to mimic the Transformation-Optics- 
type (TO-type) field transformation between two 
coordinate systems that leave the electromagnetic 
fields intact. Our primary motivation was to see if 
we could mitigate the problems encountered 
during the realization of materials prescribed by 
the TO-based designs, which are anisotropic and 
which vary over a wide parametric range. We 
have found that although the material realization 
problems are less formidable with the WEB 
approach because the materials we need are 
isotropic, the WEB approach still faces realization 
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Fig.28. Flat concave lens: (a) side view; (b) top view 
and (c) simulation result. 


Table.1. Relative permittivity values of rings of the flat 
concave lens. 


Ring 
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problems and/or matching issues. Finally, we have 
presented an alternative to the TO-type of field 
transformation approach, which is based on the 
General Scattering Matrix (GSM) concept, and 
have shown how it can help us resolve the 
material realization issue provided we modify the 
performance metric of the electromagnetic devices 
from that employed by the TO-type algorithms. 
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